function plot_gould_met5m(t1,t2)
% plot_gould_met5m(t1,t2) plots the 5-min
% basic met data in gould_met5m for the 
% period [t1:t2]. t1,t2 are optional input.
% default is t1=min(t), t2=max(t).
% 4/11/01

% get data in t1-t2
load gould_met5m
if nargin==0
   t1=min(t);t2=max(t);
end
if nargin==1
   t2=t1+1;
end
m=find(t>=t1 & t<=t2);
t=t(m);at=at(m);bp=bp(m);sst=sst(m);ws=ws(m);
wd=wd(m);sw=sw(m);lw=lw(m);rh=rh(m);
% setup figure
clf
orient tall
t0=floor(min(t));t1=ceil(max(t));
T0=floor(min([min(at) min(sst)]));
T1=ceil(max([max(at) max(sst)]));

subplot(6,1,1)
plot(t,ws,'b')
ylabel('WS')
grid
axis([t0 t1 0 30]);
title('R/V Gould 01-03 5-min Meteorological Data')

subplot(6,1,2)
plot(t,wd,'b')
ylabel('WD')
grid
axis([t0 t1 -200 200]);

subplot(6,1,3)
plot(t,at,'b',t,sst,'r')
ylabel('Ta, Ts(red)')
grid
axis([t0 t1 T0 T1]);

subplot(6,1,4)
plot(t,rh,'b')
ylabel('RH')
grid
axis([t0 t1 50 110])

subplot(6,1,5)
plot(t,bp,'b')
ylabel('BP')
grid
axis([t0 t1 925 1025]);

subplot(6,1,6)
plot(t,sw,'b',t,lw,'r')
ylabel('SW, LW(red)')
grid
axis([t0 t1 0 500]);

xlabel('2001 Yearday')
hold off